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Abstract 

This is a review of what is currently known about the gluon Sivers distribution and of what 
are the opportunities to learn more about it. Because single transverse spin asymmetries in —)■ 
TT X provide only indirect information about the gluon Sivers function through the relation with 
the quark-gluon and tri-gluon Qiu-Sterman functions, current data from hadronic collisions at 
RHIC have not yet been translated into a solid constraint on the gluon Sivers function. SIDIS 
data, including the COMPASS deuteron data, allow for a gluon Sivers contribution that is of the 
natural size expected from large N^ arguments, which is 0{1/Nf,) times the nonsinglet quark Sivers 
contribution. Several very promising processes to measure the gluon Sivers effect directly have been 
put forward, which apart from ongoing and future investigations at RHIC, would strongly favor 
experiments at AFTER@LHC and a possible future Electron-Ion Collider. Due to the inherent 
process dependence of TMDs, the gluon Sivers TMD probed in the various processes are not 
necessarily the same, but rather can be different linear combinations of two universal gluon Sivers 
functions that have different behavior under charge conjugation and that therefore satisfy different 
theoretical constraints. Eor this reason both hadronic and DIS type of collisions are essential in 
the study of the role of gluons in transversely polarized protons. 

1 The Sivers function and its definition 

The distribution of quarks and gluons in a proton (or any other spin-1/2 hadron) that is polarized 
transversely to its momentum need not be left-right symmetric w.r.t. the plane spanned by the 
momentum and spin directions. This asymmetry is called the Sivers effect [1]. It results in angular 
asymmetries of produced particles in high energy scattering processes involving a transversely 
polarized hadron. Experimental data in support for such a left-right asymmetry in the quark 
distribution was hrst obtained from semi-inclusive DIS process by the HERMES collaboration [2]. 
This review is about what is currently known about the gluonic Sivers effect distribution. The Sivers 
effect is of great interest theoretically as it is very sensitive to the color flow in the scattering process 
and to the multitude of color exchanges among initial and final states. It is the hrst quantity for 
which this has been recognized and for which color how sensitivity can be tested unambiguously [3] . 
Verihcation of its unusual properties will provide a strong test of the formalism of transverse 
momentum dependent parton distributions. It is an important quantity of nonperturbative QCD 
to consider both qualitatively and quantitatively. This review discusses these aspects for the gluon 
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Sivers distribution specifically. We first start with its proper definition. 

The number density in momentum space of a generic parton (quark, antiquark, gluon) inside a 
hadron with mass M, transverse polarization St and momentum P can be written as 

f{x, kj_-, St) = fiix, fcl) - ^ /i^(x, k'i) , (1) 

where fi{x,k‘^) is the unpolarized Transverse-Momentum Dependent (TMD) parton distribution, 
St = and P = P/\P\. The function f^{x,k‘j_) describing the distortion in the distri¬ 

bution of unpolarized partons with light-front momentum fraction x and transverse momentum 
k± due to the transverse polarization of the hadron is called the Sivers function. The notation 
used here comes from [5], but also the notation from [5l [6] is sometimes used, where 

A^ fg/ht = —2 (analogous to the quark case [7]). 

The Sivers function satisfies the following positivity bound [8]: 

^i)l ^/i(^) ^i)- (2) 

In Ref. [H] the operator definition of the gluon Sivers function, which was called Gt = [S]) was 

first given without gauge links. The definition including gauge links then appeared in Refs. mm 


^ /ir ^i) = \ /(®> St) - f{x, fc_L; -St) 


(3) 


with 




f{x,k±;ST) = 


dz- 


d2 


Ak-z 


xP+ 


(27r) 


(P, St\ 2Tr[p+^'(0)C/[o,,]F+'(z)[/[,,o]l St) 


z +=0 


(4) 


where P[a,b] is a Wilson line connecting the points a and b along a contour determined by the 
physical process and = xP^ is the fraction of parton light-front momentum. For a proper 

dehnition that is free from rapidity divergences associated with gauge links with paths (partly) 
along the light front, a redehnition involving the so-called soft factor is necessary mmmm- 
This will however not play a significant role here and can simply be considered as implicit. 


2 Sivers effect, A^, and Qiu-Sterman effect 

The Sivers effect (for both quarks and gluons) was first suggested in [T] as an explanation for the 
large left-right single transverse spin asymmetries (M^r) observed in p ^ t X [TUI fT71 fTFl nUl [2U1 
[2Tll22l[23] (and similar asymmetries in K |22j . in p |24j . and tentatively in J/^jJ [25] production). The 
Sivers effect was first studied phenomenologically in [5]. Extraction of the Sivers TMD presumes 
all-order TMD factorization, however. A factorized description of the process pp —)■ vrX only 
applies for large transverse momentum pT of the produced pions (say for GeV), where in 

fact collinear factorization is appropriate [26|, rather than TMD factorization. Although TMDs 
do appear in the phenomenological description of in [5] and subsequent studies (see m). that 
description is thus not based on a TMD factorization theorem. Rather it should be considered as an 
effective model description, now commonly referred to as the Generalized Parton Model (GPM). As 
a phenomenological approach it has proven useful in the quest to disentangle the possible underlying 
mechanisms of the spin asymmetries, see |27l [28] for more discussion, but the extracted “effective” 
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TMDs may differ from the TMDs extracted from TMD factorizing processes. What is known about 
the effective gluon Sivers TMD will be discussed below. 

In collinear factorization the single spin asymmetry (SSA) will arise at the twist-3 level [29]. In 
this description probes the (quark-gluon) Qiu-Sterman functioiii] [30] [3T] 


Tq^F{x,x) = 


M r dz dr? (P X StV 
P+ J 2-k ^ 2M 




| 2 ±l= 0 ’ 


(5) 


and its tri-gluon correlation analogues Tq\x,x) = T^\x^x) and T^g\x,x) = \x,x) [3211331 

[Ml [35] 


T^^\x,x) = 

2M5!p 
~ xiP+y 


dr] {P x StV 

2-71 2M 
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( 6 ) 


where 0^“^ = and 0^”^ = and where the light-front gauge A'^ = 0 has been considered 
for convenience. 

In |36] the quark-gluon Qiu-Sterman function has been related to the first transverse moment 
of the quark Sivers function, i.e. oc Tq^F{x,x)/M, where 


f-L(i)g 

JiT 


[x = 




(7) 


However, that relation was only established at tree level (beyond tree level the relation will be 
affected by the considered regularized definition of the Sivers TMD including its dependence on the 
soft factor). A similar tree level relation can be established in the gluon sector as well: oc 

Tg{x,x)/M (which Tq appears in this relation depends on the gauge links, see the discussion in 
section [5]) . 

Another relation has been established in m- 


fyy{x,k\) 


M ^ , , , 

~ as-^ {K ®Tq^F){x), 

fc I 


( 8 ) 


which means that the quark-gluon Qiu-Sterman function determines the large transverse momen¬ 
tum tail of the quark Sivers function. Here it should be emphasized that the function [K ® Tq^F) (x) 
consists not only of a convolution of Tq^F{x,x), but also of its derivative xdTq^F{x, x)/dx and of 
the more general Tq^F{x-,y) with y / x. It corresponds to the fact that the evolution of Tq^F is 
non-autonomous and inhomogeneous, see [38] • Note that here we have discussed the non-singlet 
contributions only, that apply to combinations like u minus d quarks, otherwise also gluonic con¬ 
tributions need to be taken into account. 

Similarly, the tail of the gluon Sivers function is determined by several Qiu-Sterman func¬ 
tions [39]. It receives contributions from the quark-gluon Qiu-Sterman functions Tq^F{x,x) and 
Tq^F{x,y) with y / X, and from the tri-gluon functions Tq^^\x,x) (which one(s) depends on the 
gauge links, see section [5]). At small x the situation simplihes: the contributions from Tg^i?(x,x) 
and Tq^F{x,y) with y 7 ^ x to the tail of the gluon Sivers function cancel each other [39]; leaving 

^Sometimes one factor of the coupling constant g is included in the definition of Qiu-Sterman functions, because 
one always encounters them multiplied by g. 
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only the tri-gluon correlators. Moreover, Tq\x,x) evolves with the same 1/x behavior at small x 
as the unpolarized gluon distribution, and is therefore not necessarily suppressed at high energies 
and small values of x, whereas Tq {x, x) lacks this 1/x enhancement [35]. As the large pr An data 
from RHIC are generally not in the small-x region of the polarized proton, except for negative xp, 
any simplifications at small x should of course first be tested for validity. 

Information from An measurements at sufficiently large pT (in order to consider a collinear 
factorization description in the hrst place) can thus in principle reflect some information on Sivers 
functions (i.e. on the tails and perhaps also on first transverse moments), but in practice other twist- 
3 contributions beside the mentioned Qiu-Sterman functions, namely chiral-odd and fragmentation 
function analogues, contribute to Aat gomiiiiaEHi. From the smallness of A^v in the midrapity 
and backward (negative xp) regions, one would generally conclude that gluonic and sea quark 
contributions to the transverse single spin asymmetry are not large, but a detailed analysis is 
required to determine precisely the size of the various contributions. 

An for TT^ production at midrapidity has been measured by the PHENIX experiment in polarized 
pp collisions at RHIC and was found to be consistent with zero, for pp values below 5 GeV at the 
permille level and for higher pp values (up to 11 GeV) at the few percent level |21[ 143 j. These data 
taken at y/s = 200 GeV probe x values only down to x ~ 0.006, where still a combination of Qiu- 
Sterman functions is expected to contribute. In [44| these vr^ data were discussed, using two models 
for the tri-gluon Qiu-Sterman functions that were constrained from A^r in H-meson production [451 
I46j . The midrapidity vr^ data are shown to mostly constrain Tq\x, x) (their N{x)) at low pp. The 
authors conclude that “Both models give tiny asymmetry due to the small partonic cross sections, so 
the form of the three-gluon correlation functions is not much constrained by the data in this region.” 
From this limited model study of both vr^ and D production one would conclude that Tq\x, x) and 
Tq\x,x) are in any case small, a permille fraction of x times the unpolarized gluon distribution. 
This should be investigated further with more general model forms that adhere to the correct 
small-x behavior and with more precise data. The experimental precision of A^v can be improved 
much further both at RHIG and especially at the AFTER@LHG experiment [39|, which would 
have a luminosity factor of 10-100 higher, if not more. Such improvement and the measurements 
of asymmetries for many different types of produced particles are required to separately constrain 
or determine the various Qiu-Sterman contributions. 

In the GPM the smallness of A^r at midrapidity puts strong constraints on the effective gluon 
Sivers function. As explained, this gluon Sivers function captures the combined effect from several 
Qiu-Sterman contributions and may thus differ from the actual gluon Sivers function obtained 
from TMD-factorizing processes. In a recent GPM analysis HZ], which is an updated analysis 
of Ref. [38], the best fits to the PHENIX midrapidity vr^ Ajv data indeed correspond to a small 
effective gluon Sivers function w.r.t. its theoretical bound determined by the unpolarized gluon. 
Eor example, for x < 0.1 it is at most only a few percent of the bound. However, the maximally 
allowed effective gluon Sivers function is still sizable though. Its first transverse moment is still 
found to be around 30% of the up quark Sivers function in the region 0.06 < x < 0.3, which is 
consistent with findings from semi-inclusive DIS for the actual gluon Sivers function and also with 
theoretical expectations, as discussed in the next section. In addition, it should be mentioned that 
this GPM analysis assumes a Gaussian k± dependence, which does not correspond to the correct 
power-law tail of the Sivers function, Eq. ([8|), nor of the unpolarized gluon distribution. Given all 
the caveats that come with these results, one should be careful to draw a definite conclusion about 
the size of the actual gluon Sivers TMD from An data. 

One should also specify clearly what one calls a small gluon Sivers function. It will depend 
on what one compares to, i.e. whether that is to the unpolarized gluon that grows very rapidly at 
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small X, or to the up or down quark Sivers function for not too small x. At small x it becomes very 
important whether one discusses the f or d type contribution (see section ED, which is an issue not 
addressed in the GPM studies of 

3 Sivers asymmetry in SIDIS 

The Sivers effect leads to a — asymmetry in semi-inclusive DIS (SIDIS) [4], which has been 

observed in experiments using a proton target by HERMES [21 [50] and COMPASS [51], and using a 
^He target by Jefferson Lab Hall A [52l|53]. The data follow to quite a good extent the expectations 
of a valence quark picture in the target and of favored fragmentation. In the proton case, the Tr"*" 
thus shows the largest asymmetry (for large z values, the asymmetry is around 4-5%, or even 
somewhat larger when a lower cut of > 4 GeV^ is implemented instead of > 1 GeV^ [50]). 
The Tr~ asymmetries are smaller and still compatible with zero. The asymmetries are similar 
to the TT^ asymmetries, but with larger errors. Sivers asymmetries on a deuteron target [M] are all 
consistent with zero. Fits to all these HERMES and COMPASS data, including the deuteron data 
using isospin symmetry, indicate that the Sivers function (/f^) for u quarks in a proton is negative 
and for a d quark in a proton is positive and approximately equal in absolute value [55]. This fits 
the expectations from the limit of a large number of colors Nc |56[ 157] : 

kl) = kl) + 0{1/N,). (9) 

The flavor singlet combination of u and d is of the same order as the gluon contribution in Nc 
counting [58]. The latter is thus 1/Nc suppressed w.r.t the flavor non-singlet quark Sivers effect at 
not too small x (x ~ 1/A"c) [59] . 

Within the current accuracy, the SIDIS data do not require any sea quark or gluon contributions, 
which among other considerations (see section [6]) led Brodsky and Gardner to conclude that the 
gluon Sivers function is small or even zero (“absence of gluon orbital angular momentum”) [60] . The 
SIDIS data from HERMES, COMPASS, and Jefferson Lab Hall A are of course at rather modest 
and not too small x values, i.e. in the valence region. One cannot yet draw any conclusions about 
the gluon Sivers function at higher and smaller values of x. Moreover, the data certainly still 
allow for gluon Sivers contributions of the order of 1/Nc times the valence quark Sivers functions. 
This is evident from the fits by Anselmino et al. [55], where the first transverse moment of the u 
and d Sivers functions have error bands that are at least around 30% of the central values. 

Note that the SSA in the “inclusive” process ep —)• hX, where the back-scattered lepton is not 
observed [61(162]. does not allow for an interpretation in terms of TMDs, as the data are dominated 
by ~ 0. Even for large pT the appropriate factorization would be collinear factorization and the 
Sivers type of asymmetry would probe the Qiu-Sterman functions instead [63], which as discussed 
above have some relation to the Sivers TMDs, but only via the tail or possibly via the first transverse 
moment. The asymmetries for > 1 GeV are found to be at the level of 5-10% for positive hadrons. 
Fits will need to make clear how much room there is for a gluon Qiu-Sterman effect. Given the 
fact that the gluon Qiu-Sterman function does not enter at leading order in as in this process, this 
room may be considerable. 

4 Sivers asymmetry in other processes 

Several other pp scattering processes to access the gluon Sivers function have been suggested over 
the past years: p"^ p —)• jet jet X [ 6 l ] , p'^ p ^ D X [65] (M] [M] , p''' p —)• 7 X [ 66 ] , p"^ p —^ 7 jet X [HQ [70] , 
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p —)■ 7 *X —)• p~ X [ 66 ], p^p —)• p^p —)■ vrjetX [^, p"^p Vc/b^ [SS]. Several of 

these processes are like ^at in (high-pT) pion production, which means that they deal with twist-3 
collinear factorization and only provide indirect or limited information about the gluon Sivers TMD. 
Several other processes run into the problem of TMD factorization breaking contributions [72| and 
hence are not safe. In principle they do probe TMDs but as a result of TMD factorization breaking 
contributions, conclusions about the gluon Sivers function from their measurements cannot be 
drawn safely. This applies for instance to the process —)• jet jet X (measured at RHIC to 
be small at the few percent level |73jl. which moreover suffers from cancellations between u and 
d contributions and between the effects of initial and final state interactions [711|75l|76|. TMD 
factorization breaking would also apply to open heavy quark production: p^p ^ QQX, such as 
p^p —)• D^D^X, cf. e.g. |77j . Whether the problem also applies to double heavy quarkonium 
production remains to be seen, because in practice the color singlet contributions may give the 
dominant contribution in that case. Among the hadronic collisions the processes having one or 
two color singlets in the final state would in any case be safest. One very promising example is 
p^p —>• 7 jet A [TO], where it depends on the rapidity of the photon and the jet, i.e. on the x 
fraction of the parton in polarized proton, whether the gluon Sivers function dominates over the 
quark one or vice versa. Another very promising example isp^p—J/V'qA, which is predominantly 
initiated by gluon-gluon scattering which is an order in as higher than the gluon contribution in 
p^p —)• 7 jet A) and for which the color singlet contribution dominates over the color octet one to a 
large extent [781179]. The same applies top^p^ J /tp J/il^ X (see the contribution by Lansberg and 
Shao in this special issue). AFTER@LHC would be very well suited for studying these processes. 

SSA experiments could be done at AFTER@LHC where the beam of protons or lead ions of the 
LHC would collide with a fixed target that is transversely polarized. Such pp^ and Pbp^ collisions 
would have a center-of-mass energy of 115 and 72 GeV, respectively, and have high luminosity 

and good coverage in the rapidity region of the transversely polarized target (mid and large Xp) [39]. 
Polarized Drell-Yan and prompt photon production studies could be done to measure the quark 
Sivers function very precisely, perhaps to the level that the gluon Sivers function becomes relevant, 
despite the large values of x in the polarized target. As mentioned yjet and J/ip'j production could 
be used to study the gluon Sivers effect directly, where the former would need specific selection of the 
rapidities. In addition, the comparison of Pbp^ —)■ 7 jet A and pp^ —?> 7 jet A would give a further 
handle on determining the relative sizes of quark and gluon Sivers functions. Other processes, such 
as D-meson or J/pi production, would allow a similar study of Qiu-Sterman functions, including 
the tri-gluon ones, which are of course interesting in their own right. See [28] for a more detailed 
and quantitative study of twist-3 transverse single-spin asymmetries in proton-proton collisions 
at the AETER@LHC experiment. All these possibilities offer a very interesting complementary 
opportunity or even a competitive alternative to the other existing high-energy particle physics 
spin projects aiming at studying the role of gluons in transversely polarized protons. 

In electron-proton scattering one of the most promising processes to directly probe the gluon 
Sivers function is open charm production, ep^ — e' ccX, which could ideally be studied with an 
Electron-Ion Collider (EIC). By selecting the charm (or bottom) quark, one effectively eliminates 
the subprocesses 'y*q —)• qg and 7*5 —)• qg and becomes essentiall}!! sensitive to 7*5 —cc, and thus 
to the gluon Sivers functior0. Here the transverse momenta of the heavy quarks are considered to 

^Single transverse spin asymmetries in jet production measured at RHIC |67II68| at forward rapidities (the valence 
region) show very small asymmetries, which is probably due to a cancellation among u and d quark contributions 

m- 

®This assumes that intrinsic charm contributions are suppressed by selecting sufficiently small x values. 

^A similar argument is used in the study of high-pT hadron pairs in muon-deuteron and muon-proton scattering |80l 
I81| . where photon-gluon fusion is expected to dominate. The relevant asymmetry is found to be —0.14± 
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be almost back-to-back. There is no problem with TMD factorization breaking contributions of the 
type discussed in m, but that does not mean the process is as straightfoward as SIDIS. Even in the 
case where one considers charm jets, one has to include a description of the transverse momentum 
distribution inside such a jet. It may be easier to consider measurements (for a study of the 

twist-3 SSA in large pT D meson production in SIDIS, i.e. ep^ —)• e' DX, see [SSlES]). In either 
case one deals with 3 TMDs. Such processes involve a different soft factor (in this case a vacuum 
correlator of 6 Wilson lines) compared to processes involving 2 TMDs as in SIDIS, affecting the 
predictability. This has been discussed at the one-loop level in [83]. The SSA in —)• e' X 

has been studied for some models of the gluon Sivers function in [ 8 l|, cf. section 2.3.1. This may be 
the ‘smoking gun’ process for the gluonic Sivers effect at an EIC. It should be mentioned though 
that it actually probes a different gluon Sivers TMD than the hadronic processes discussed above. 
This is discussed in the next section. It shows that hadronic processes are complementary to DIS 
processes. 

For completeness we mention that when comparing extractions of the gluon Sivers TMD from 
different processes, one not only has to take care of the process dependence, but also of the dif¬ 
ferent energy scales. Under TMD evolution from one scale to another, the transverse momentum 
distribution changes. For details we refer to [85l [ 86 l [571ES 15^ 190] . 


5 Process dependence of the gluon Sivers function 


Once a set of processes that in principle allow to probe the gluon Sivers TMD has been obtained, 
one still has to take into account the fact that such a TMD is process dependent. For quarks the 
famous overall sign change between the Sivers TMD probed in SIDIS and the one probed in Drell- 
Yan is expected [HU |92l[3l[93| and is currently under experimental investigation. For gluons the 
situation is more complicated as each gluon TMD depends on two gauge links (in the fundamental 
representation), so there are more possibilities [9ll[95l[TT|. The gauge link structure of the gluon 
distributions ui ep ^ e' X differs from the one in for instance pp —^ 7 jet A (cf. [96] for the 

comparison at small x). Clearly, this will complicate the analysis of gluon Sivers effect which will 
involve more than one gluon Sivers function. In HI] it was demonstrated that any gluon Sivers 
function can be expressed in terms of two “universal” gluon Sivers functions, 


f^9[U] 

JIT 




C=1 




( 10 ) 


where the coefficients are calculable for each partonic subprocess. The first transverse moments 
of the two distinct gluon Sivers functions are related (at least at tree level) to the two distinct tri¬ 
gluon Qiu-Sterman functions . Therefore, we will refer to the universal gluon Sivers functions 

as and ■ The two functions have different behavior under charge conjugation, just 

like is a matrix element of a C-even operator and of a C-odd operator. 

The process ep^ —)• e' X is dominated by just one partonic subprocess 75 ^ qq and thus 

probes the gluon Sivers function with two future-pointing (-I-) links |97j . which is jllj . The 

process p^ p —^ 7 jet A probes the subprocesses qg —>• jq and qq —)• 75 . If one selects kinematics 
such that one probes small x values in the polarized proton, such that qg —)• jq dominates, then 

0.15(stat.) ± 0.06(syst.) at {xq) = 0.13 for the deuteron and —0.26 ± 0.09(stat.) ± 0.08(syst.) at (xq) = 0.15 for the 
proton. For the interpretation of the data in terms of the gluon Sivers effect, and the pr of each hadron need to 
be sufficiently large to trust factorization. 
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this process accesses the gluon Sivers with a future and a past-pointing link, which corresponds to 
. The theoretical expectations are different for these two cases. 


6 Theoretical constraints on Sivers function 


Constraints on the unintegrated gluon Sivers TMD f^{x, from fits have to take into account 
that it is theoretically possible that both the quark and the gluon Sivers TMD can have nodes 
in X and/or k± [98l [99]. The possibility of a node in x is supported by the observation [35| that 
the splitting function for Tq is negative at small x, in analogy to the Ag case. Fits to SIDIS 
data (studied with a rather restrictive parameterization and in a restricted kinematic range) do 
not appear to require a node |100j . but that does not exclude this possibility. Especially when 
comparing data from different kinematic regions and different processes, this option should be kept 
in mind. Nodes can of course have a large effect on integrals of Sivers functions, such as the first 
transverse moment (© and its hrst Mellin moment (for parton a) 

{k±a) = -M{St xP)jdx (11) 


which is the average transverse momentum inside a transversely polarized target. The notation 
(fcj_a) comes from [lOlj . This quantity is related to the Sivers shift |lf)2] . the average transverse 
momentum shift orthogonal to the transverse spin direction, which is normalized to the zeroth 
transverse moment of the unpolarized TMD /^^^(x) = f d?kx /i(x, k\): 


{k\{x))uT^ 


= M 




( 12 ) 


Here only the y-component perpendicular to the transverse spin direction x is nonzero and therefore 
considered. Note that the Sivers shift depends in principle on the renormalization scale /r, a rapidity 
variable C; but also on the path of the gauge link (see below). Furthermore, there is the theoretical 
issue whether the transverse moments converge. For that reason a regularized version using Bessel 
moments has been suggested in |lf)2j . which for quarks has been evaluated on the lattice in [103] . 
The lattice calculation of the Sivers shift for u — d are consistent with negative u and positive 
d Sivers functions in SIDIS, which are expectations that follow from general arguments on final 
state interactions nnn, from a model-dependent relation to GPDs [10511106] . see below, and from 

fits [ToTirms] . 

As said above, f^^^'^{x) has a (process-dependent!) relation to the Qiu-Sterman function 
Tg^p{x, x), which so far only has been established at tree level and for quarks. In addition, Burkardt 
has suggested a model-dependent relation between the integrated quantity (x) and an integral 
over the GPD E(x,^, A^) at zero skewness .^ = 0 (and only for quarks) [105] : 

fi~T^\x)(x J d‘^b_iT{b±)-^£{x,b\), (13) 

for a nucleon polarized in the transverse x direction. Here T(x, b']_) = f J 2 ^t E(x, 0, — A^) 

andZ(b±) is called the leasing function. This relation has been obtained in models |105[ 1106] 1^ 1109] . 
It allows to relate (kpg) to the anomalous magnetic moment Kq associated with the quark q: 

























albeit in a model-dependent and (due to the different integrals involved) only qualitative way. This 
relation does confirm the expectations for the relative signs between the u and d Sivers functions, 
and has been used in Ref. m to fit SIDIS data for the Sivers effect with the integral constrained 
by the anomalous magnetic moments. Interestingly, this led to a new estimation of the quark total 
angular momentum which turned out to be in agreement with most common GPD extractions m 
11121 11131 11141 1115] . The relation between (k^g) and Kg is also at the heart of the argument 
by Brodsky and Gardner of why a gluon Sivers function is expected to be small. Using that 
kZ = 2Kp + Kn = 1.673 and = 2k„ + Kp = —2.033, one sees the opposite signs reflected, but 
since |(Kp + Kn)/2\ = 0.06 <C Kp/n ~ 1-8 — 1.9 this suggests that there is little room for gluon 
contributions |60j . If (kS + /2 is taken as a measure for Kg, the latter is about 10% smaller than 

Kg. This would suggest that k^ is of order l/N^ rather than 1/Nc, which in turn would suggest a 
similar additional 1/Nc suppression for the gluon Sivers function. Glearly there are various (strong!) 
assumptions going into this type of argument, such that the conclusion can certainly not be taken 
at face value. Apart from the assumptions on the relation to the gluon Sivers function, it is not 
clear that one can use very low-energy quantities to deduce something about the size of the gluon 
contributions at energies around or above 1 GeV to begin with. 

Burkardt derived a further constraint on the fully integrated quantity {k±a), nowadays referred 
to as the Burkardt sum rule (BSR) [101] . stating that the total transverse momentum of all partons 
in a transversely polarized proton must vanish m, 

(fcx)= E if^±a) = 0. (15) 

a.=q,q,g 

In terms of the Sivers function, the BSR takes the form |59j 

f dx = 0. (16) 

a=q,q,g '' 

Its validity has been checked explicitly in a diquark spectator model in Ref. 
data from Ref. [108] at the scale = 2.4 GeV^ almost saturate the BSR 
d quark contributions alone: 

{k±u) = 96^®° MeV, (k^d) = -1131^? MeV. 

The contributions of the sea quarks are all small and together allow the 
gluon contribution: 

- 10 < {k_ig) < 48 MeV. (18) 

This means there is certainly still room for a 30% contribution from gluons w.r.t. the valence quarks. 
Of course, it should be emphasized that these values were obtained under assumptions on the k± 
dependence, the absence of nodes, and extrapolations outside the kinematic region accessed by the 
SIDIS experiments. 

The derivation of the BSR by Burkardt [101] considers gauge links as appear in SIDIS and 
involves a gluon correlator containing the antisymmetric fate structure constant of SU{3). As 
shown in [941195] . there is also a gluon correlator with the symmetric date structure constant. As a 
consequence, inclusion of gauge links in the operator definition of TMD distributions gives rise to 
two distinct gluon Sivers functions, and (corresponding to the (Ac) label used in [TT]'). 

However, the BSR essentially expresses transverse momentum conservation. Since the momentum 
operator in QCD is C-even, only the gluon Sivers function which is associated with a C-even 


[TTT] . The fits to SIDIS 
already with the u and 

(17) 

following range for the 
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operator is constrained by the BSR. The gluon Sivers function , which is associated with a 

C-odd operator, is not expected to satisfy a BSR where quark and gluon contributions cancel each 
other. Judging from the small x behavior of the Tq\x,x) expected from its evolution equation, 
the integral of over x may even not converge. 

7 Conclusions 

In summary, no hard constraints on the size of the gluon Sivers function exist apart from the 
positivity bound, although the theoretical expectation from large considerations (expected to 
hold approximately for not too small x) favors a 30% gluon to quark Sivers ratio, which is still 
completely allowed by all SIDIS data, including the COMPASS deuteron data. It may turn out 
that the ratio is smaller, but much smaller than 10% may in turn be considered unnaturally small. 
Strictly speaking, no direct conclusion about the size of the gluon Sivers function can be drawn 
from An data. Like p —)• tt X, many other processes suggested in the literature to probe the 
gluon Sivers function actually deal with collinear factorization and as such they are sensitive to 
complicated linear combinations of quark-gluon and tri-gluon Qiu-Sterman functions (and chiral- 
odd and fragmentation function versions of them) rather than to Sivers functions directly. Inferring 
constraints on the gluon Sivers function, even on its large transverse momentum tail must therefore 
be done with much care. In the Generalized Parton Model description of Atv at midrapidity, the 
effective gluon Sivers function is currently still allowed to be 30% of the up quark Sivers function, 
despite the smallness of the asymmetry. Other suggested processes that in principle probe TMDs 
may suffer from TMD factorization breaking contributions and any results on the gluon Sivers 
function from measurements of such processes cannot be trusted. This applies for instance to the 
process ^ jet jet A. 

The most promising processes that directly give access to the gluon Sivers effect are p^p —)■ 
7 jet A, p^p —)• J/V'yA and ep^ —e'ccA. The hrst process can be studied at RHIC and at a po¬ 
larized hxed-target experiment at LHC (AFTEROLHC), the second process also at AFTER@LHC, 
and the third process at a possible future Electron-Ion Collider. Due to the inherent process de¬ 
pendence of TMDs, the gluon Sivers TMD probed is in principle different in these processes. They 
can be expressed in terms of two universal gluon Sivers functions that appear in different linear 
combinations in different processes. Extracting and comparing these universal functions is very 
interesting from a theoretical point of view. The fact that a difference can exist is a consequence of 
the non-Abelian nature of QCD. Both functions satisfy different theoretical constraints. Although 
TMD factorization is expected to hold for these processes, that has not been demonstrated yet to 
all orders. Apart from the process dependence, there is also the issue of modified soft factors to 
contend with still. Nevertheless, as far as experimentally demonstrating and measuring a gluon 
Sivers effect in transversely polarized protons, several complementary future possibilities exist, in 
which AFTER@LHC can play a very important role. 
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